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Abstract 

We study the thermodynamics of the Bose-condensed atomic hydrogen con- 
fined in the loffe-Pritchard potential. Such a trapping potential, that models 
the magnetic trap used in recent experiments with hydrogen, is anharmonic 
and strongly anisotropic. We calculate the ground-state properties, the con- 
densed and non-condensed fraction and the Bose-Einstein transition temper- 
ature. The thermodynamics of the system is strongly affected by the anhar- 
monicity of this external trap. Finally, we consider the possibility to detect 
Josephson-like currents by creating a double-well barrier with a laser beam. 
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Few years ago, Bose-Einstein condensation (BEC) has been experimentally observed in 
clouds of trapped alkali-metal atoms [1] . Recently, BEC has been also achieved with atomic 
hydrogen confined in a loffe-Pritchard trap [2]. That is an important result because hy- 
drogen properties, like interatomic potentials and spin relaxation rates, are well understood 
theoretically. As stressed by Killian et al. [2] , the s-wave scattering length of the hydrogen 
is very low and, compared with other atomic species, the condensate density is high, even 
for small condensate fractions. Moreover, due to hydrogen's small mass, the BEC transition 
occurs at higher temperatures than in those of alkali atoms. 

In this paper we calculate the thermodynamical properties of the trapped hydrogen gas 
by using the quasi-classical Hartree-Fock approximation [3]. This approach is justified by 
the very large number of atoms (about lO^*^) in the trap and by the relatively high temper- 
atures involved (order of /xK). Due to the anharmonic external trap, the analytical results 
for BEC thermodynamics obtained by Stringari et. al [3] cannot be used for quantitative 
predictions. Our detailed theoretical study of the hydrogen thermodynamics can give useful 
informations for future experiments with a better optical resolution. In the last part of the 
paper we discuss the criteria for macroscopic quantum tunneling and macroscopic quantum 
self-trapping by using a laser beam to create a double- well potential. 

In the experiment reported in Ref. [2], the axially symmetric magnetic trap is modelled 
by the loffe-Pritchard potential 

C/(p,^) = VM' + (/^^' + 7)'-7, (1) 

where p and z are cylindrical coordinates and the parameters a, (3 and 7 can be calculated 
from the magnetic coil geometry. In particular, for small displacements, the radial oscillation 
frequency is uip = a/y/rnj = 27t x 3.90 kHz, the axial frequency is uj^ = ^2/3/m = 2tt x 10.2 
Hz and -f/ke = 35 //K [2]. 

A dilute gas of N identical hydrogen atoms is theoretically described by a bosonic field 
operator ■?/'(r, t). Following a standard treatment, one can separate out the condensate part 
with the Bogoliubov prescription tp{r,t) — $(r) -|- ^(r, i), where <l>(r) — {ip{r,t)) is the 
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macroscopic wavefunction (order parameter) of the condensate, normalized to the number 
A'o of condensed atoms, and (...) is the mean value in the grand-canonical ensemble. Then, 
in the mean-field approximation, the order parameter $(r) satisfies the following finite- 
temperature Gross-Pitaevskii equation [4,5] 



-— + U(r) + gno(r) + 2gnT(r) 



$(r) = //$(r) , (2) 



where m is the mass of the atomic hydrogen, /j, the chemical potential, and g — Anh^a/m 
is the scattering amphtude with a the s-wave scattering length (a — 0.0648 nm). The 
condensate density is no(r) = |$(i')P and nT(r) =< 0"'"(r, t)0(r, t) > is the thermal density 
of non-condensed particles, normalized to Nt — N — Nq. Notice that we ignore the T = 
quantum depletion. In fact, the corrections to the mean-field results are fixed by the gas 
parameter n(0)a^, where the total density n(r) is evaluated at the center of the trap [6]. 
We shall show that this parameter is very small also when there are billions of atoms in the 
trap. 

The thermal density can be calculated through the quasi-classical Hartree-Fock approx- 
imation 

nrir) = 1 ^3/2 (e-(^W+^^"«-/^)AB^) , (3) 

where 510,(2;) = Efeli-^*^/^" and A = /mkBTY^'^ is the thermal length. Thus, the 

thermal particles behave as "non-interacting" bosons moving in the self-consistent effective 
potential U{r) -\- 2gn{r), where the term 2gn{v) is the mean-field generated by interactions 
with other atoms [3] . Note that the quasi-classical approximation is accurate for the exper- 
imental hydrogen cloud, which has = lO^'^ atoms. Also the Hartree-Fock approximation 
is valid in this context because /c^T is much larger than the chemical potential /x. 

In many papers, the thermodynamics of Bose-condensed dilute gases has been studied by 
solving the Eq. (2) and (3) with a self-consistent iterative procedure (see [6] and references 
therein). We shall adopt this method that, in our case, describes correctly both the high- 
and low-temperature regimes due to the negligible effect of collective excitations in the 



thermal density. It is important to observe that mean-field predictions for trapped alkali- 
metal atoms are in good agreement with recent path- integral Monte Carlo calculations [7]. 
Moreover, in the experiments with atomic hydrogen, the gas is dilute {na^ <^ 1) but also 
strongly interacting because Na/an ^ 1, where an — {h/muHY^'^ with ujh — {u'^pUzY^^ . It 
follows that the kinetic term of the Eq. (2) can be safely neglected, as we have verified by 
numerically solving Eq. (2) at zero temperature, and one gets the Thomas-Fermi condensate 
density 

no(r,t) = ^[/x-C/(r)-2^nr(r)] , (4) 

in the region where ji > U{y) -\- 2gnT{v), and no(r) = outside. In practice, we study the 
BEC thermodynamics by solving self-consistently the Eq. (3) and (4). 

First, let us consider a noninteracting gas {g — 0). In such a case, the formula (3) of 
thermal density simplifies and one can directly obtain the BEC transition temperature 
as a function of the total number of atoms by numerically solving the equation 

In Fig. 1 we compare the BEC transition temperature of the loffe-Pritchard potential 
with the analytic formula T°'^ = O.dAtiWHN^^^ /ks: that is exact in the thermodynamic 
limit for the harmonic potential U = (m/2)(a;^p^ -|- cu^^;^). Fig. 1 shows that, for a large 
number of atoms, T^'^ exceeds T^. In particular, for = 2 • 10^° the relative difference 
is more than 41%. In the loffe-Pritchard trap the BEC transition temperature is strongly 
affected by the fact that along the cylindrical radius the quadratic behavior becomes almost 
linear at large distances (for a discussion of potentials with exact power laws see [8]). It 
means that the density of states is higher than in the harmonic case and consequently T° is 
suppressed. The critical temperature in this loffe-Pritchard trap can be well represented by 
the law T° = bN^/"^, where b = 5.47 • lO'^ //K and rj = 3.51. 

The role of the interatomic interaction on the transition temperature Tc is very small. 
In Fig. 1 one observes that the repulsive interaction reduces T,. both in harmonic and loffe- 
Pritchard traps. The shift of the transition temperature STc is in agreement with the law 
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5Tc/Tc — — 1.3(a/aij)A^-^/^, predicted by Stringari et al. [3] for the harmonic trap. These 
corrections are of the order of 0.1 jjK. 

In Fig. 2 we show the condensate density no(r) and the thermal density nT(r) at T = 45 
and 47 jiK for N — 2.2- 10^^ atoms. The thermal density shows a depletion near the origin, 
that is due to the presence of the condensate fraction in that region. Note that the thermal 
cloud fills a very large spatial region compared to the condensate one. 

In Fig. 3 we present various properties of the condensate as a function of temperature 
for N — 2.2 ■ 10^° atoms. We consider the number of atoms and temperatures achieved in 
the MIT experiment [2] . The central density and the size of the condensate are particularly 
interesting because can be easily detected experimentally. Our results are compatible with 
the experimental data: at T 45 ± 5 jjK with N — 2.2 ■ 10^^ the estimated condensate 
fraction is 5%, the peak condensate density 4.8 ± 1.1 • 10^^ cm~^, the condensate diameter 
is d = 15 /im and its length I — 5 mm [2]. We note that, as previously anticipated, the gas 
parameter n{0)a^ is always less than 10"^. 

The condensed fraction, the energy per particle of the condensate and the chemical 
potential are shown in Table 1. We calculate also the two-body decay rate T — cf (Prnl{r), 
where c = 1.1 • 10~^^ cm^/s [3]. The decay time can be estimated as T1/2 = 7/(2cno(0)) [3]. 
The experimentally measured life-time (about 5 seconds with Nq/N = 0.05) is larger than 
Ti/2- As stressed in [2], the thermal gas continually replenishes the condensate as atoms are 
lost through two-body collisions. Thus the apparent hfe-time of the condensate is increased 
by this replenishment. 

For the sake of completeness, in Fig. 4 we plot the condensate fraction Nq/N as a, function 
of temperature for N — 10^, 10^ and 10^° atoms. Note that Hijmans, Kagan, Shlyapnikov 
and Walraven [9] have shown that, due to the balance between the thermalization rate and 
the two-body spin-relaxation decay rate, the maximum achievable equilibrium condensate 
fraction for hydrogen cannot be very large. Moreover, a recent theoretical paper [10] sug- 
gests the possibility of interesting nonequilibrium effects like the short-time formation of 
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quasicondensate droplets. However, we belive that the equihbrium results we have shown 
may provide a useful guide for future experiments with a better optical resolution. 

A very interesting issue is the possibility to detect macroscopic quantum tunneling 
(MQT) and Josephson-like oscillations of the Bose condensate in double- well traps [11,12]. 
Recently, we have shown that with ^^Na atoms in harmonic trap one sees only the macro- 
scopic quantum self-trapping (MQST) of the condensate. To get outside the MQST regime 
it is necessary to strongly reduce the scattering length [12]. Due to its very low scattering 
length, atomic hydrogen is a good candidate for MQT in double-well traps. 

It is easy to create a double-well trap for a cigar shaped condensate by using a laser 
beam [13]. The effect of a laser beam on atoms can be modelled by the following potential 



where the potential barrier Uq is proportional to the total power of the laser beam perpen- 
dicular to the long axis of the condensate, and a is the beam radius [13]. 

Smerzi et al. [11] have shown that the time-dependent behavior of the condensate in a 
double- well potential depends on the parameter A = 4£^j„(/| A£'°|: Ei^t is the interaction 
energy of the condensate and AE^ is the kinetic-|-potential energy splitting between the 
ground state and the quasi-degenerate odd first excited state of the GP equation [12,13]. 
Let z — {Ni — N2)/N be the fractional population imbalance of the condensate in the two 
wells. For a fixed A (A > 2) there exists a critical Zc — 2-\/A — 1/A such that ior < z < Zc 
there are Josephson-like oscillations of the condensate. For z <^ Zc these oscillations are 
harmonic with period r = to/vT+~/V, where Tq = 2'!rh/\AE^\. Instead for < < 1 
there is MQST of the condensate: even if the populations of the two wells are initially set 
in an asymmetric state {z ^ 0) they maintain the original population imbalance without 
transferring particles through the barrier [11,12]. 

By solving the zero-temperature GP equation (2) with the steepest-descent method in 
the double- well trap given by Eq. (1) and Eq. (6), we calculate the parameter A for the 
condensate of hydrogen atoms in the tunneling region. In Table 2 it is shown that, for a 




(6) 
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small number of condensed atoms (about 10^), the MQT is possible with a relatively large 
fractional population imbalance z. For this number of condensed atoms, the period r of 
oscillation is smaller than the life-time ri/2 of the condensate, that is about 50 seconds. 
For a larger number of atoms, the parameter A quickly grows and the condensate remains 
self-trapped. Note that at nonzero temperature, BEC depletion and thermal fluctuations 
will slightly modify the parameters of the tunneling and will damp the coherent oscillations. 
Nevertheless, as shown by Zapata, Sols and Leggett [14], the effect of damping is neghgible 
for temperatures lower than about X^hion- 

In this paper we have studied the thermodynamics of the atomic hydrogen in a loffe- 
Pritchard trap. We have found that the BEC transition temperature is strongly dependent 
on the anharmonicity of the trap while the effect of repulsive interaction is small. We have 
calculated the properties of the Rose condensate at various temperatures. Our results are 
consistent with available experimental data, which still have a large error. The present 
paper should be very useful for future experiments because it gives detailed information on 
the Bose condensate at temperatures not achieved so far. Finally, we have considered the 
inclusion of a laser beam to produce Josephson-like currents in the resulting double-well 
trap. Our calculations suggest that with a small number of condensed atoms (about 10^) 
the macroscopic quantum tunneling can be observed. 
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T(/xK) 


No/N 


E/No 


/i 


r 




44.0 


0.21 


4.64 


6.59 


2.53 


0.37 


44.5 


0.18 


4.39 


6.25 


2.08 


0.38 


45.0 


0.16 


4.13 


5.89 


1.65 


0.41 


45.5 


0.13 


3.82 


5.46 


1.27 


0.44 


46.0 


0.10 


3.45 


4.95 


0.88 


0.49 


46.5 


0.07 


3.04 


4.37 


0.56 


0.56 


47.0 


0.04 


2.52 


3.64 


0.29 


0.69 


47.5 


0.02 


1.89 


2.75 


0.10 


0.93 



TAB. 1. Condensate properties, for TV = 2.2-10^'^ interacting hydrogen atoms in the loffe-Pritchard 
trap, at various temperatures: condensed fraction Nq/N, energy per particle of the condensate 
E/Nq and chemical potential fi in units of W^hcOz (where lO^^w^/Zcs = 0.49 /xK), two-body decay 
rate F in units of 10^° s~^ and decay time ri/2 in seconds. 



N 


Uo 


Eint 




A 




ZcN 


T 


10000 


8 


1.371 


0.222 


24.69 


0.39 


3900 


0.09 


10000 


10 


1.383 


0.164 


33.64 


0.34 


3400 


0.10 


10000 


12 


1.397 


0.126 


44.38 


0.30 


3000 


0.12 


50000 


16 


4.345 


0.095 


183.03 


0.15 


7500 


0.08 


50000 


18 


4.361 


0.072 


242.95 


0.13 


6500 


0.09 


50000 


20 


4.378 


0.055 


317.80 


0.11 


5500 


0.10 



TAB. 2. Parameters of the MQT for A'^ = 10^ and N = 5 ■ 10^ condensed hydrogen atoms. 
loffe-Pritchard external potential and different values of the double- well potential barrier Uq. The 
laser-beam radius is a = 5 /xm. The energies are in units of HoOz and the oscillation period r in 
seconds. 
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FIGURE CAPTIONS 



Fig. 1. BEC transition temperature Tc versus number N of hydrogen atoms. Compar- 
ison between harmonic and loffe-Pritchard trapping potential. 

Fig. 2. Condensate density profile no and thermal density profile ut, for a cloud of 
= 2.2 • 10^° interacting hydrogen atoms in the loffe-Pritchard potential. Densities are in 
units of al- Notice the very different abscissae scales. 

Fig. 3. Observables of the Bose condensate as a function of temperature T.N — 2.2-10^° 
interacting hydrogen atoms in the loffe-Pritchard trap. 

Fig. 4. Condensate fraction Nq/N as a function of temperature T. N interacting 
hydrogen atoms in the lofi'e-Pritchard trap. 
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- loffe-Pritchard potential (non interacting) 

- Harmonic potential (non interacting) 
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